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We exploit classical results on the Stieltjes moment prob- 
lem to obtain completely explicit necessary and sufficient con- 
ditions for the photon number distribution of a radiation field 
mode to be classical. These conditions are given in two forms - 
respectively local and global in the individual photon number 
probabilities. Equivalence of the two approaches is demon- 
strated. Detailed quantitative statements on oscillations in 
the photon number probabilities are also presented. 



I. INTRODUCTION 

The study of nonclassical aspects of radiation from a 
variety of viewpoints is of considerable and continuing 
interest. On the conceptual side one can distinguish be- 
tween phase sensitive signatures of nonclassicality and 
phase insensitive ones. The most familiar of the for- 
mer is quadrature squeezing |l , while in the latter we 
have amplitude squeezing or subpoissonian photon statis- 
tics [^1^ . Higher order squeezing ||] as well as conditions 
on the factorial moments ^ of the photon number dis- 
tribution (PND) have also been presented as sufficient 
conditions for nonclassicality. In some cases, the single 
mode treatments have been generalised to the two or gen- 
eral multimode situations ItI . 

From another point of view, a finer classification of 
states of quantised radiation according to increasing non- 
classicality has been given, stressing the connection to the 
specific classes of observables being measured Q. The 
case of Gaussian- Wigner distributions for a single mode 
field has been analysed in detail to give a concrete illus- 
tration of these ideas |^ . 

In the present paper we give a complete treatment of 
phase insensitive measurements on the single mode quan- 
tized radiation field, and develop necessary and sufficient 
conditions for nonclassicality of the field. Our focus is 
on the PND, for any given state, which contains all in- 



formation concerning all phase insensitive measurements. 
There are several previous partial results in this direction. 
The best known is the Mandel Q-parameter criterion for 
distinguishing globally between sub and super poissonian 
statistics. This is global in the sense that the evaluation 
or measurement of the parameter Q in any given case re- 
quires knowledge of the probability for finding n photons 
in the state, for all n. We may also mention the impor- 
tant generalisation of the Mandel Q parameter achieved 
by Agarwal and Tara Q. Our principal motivation is 
to look for local conditions on the PND, involving only 
a finite number of photon probabilities, which will serve 
as signatures to distinguish between classical and non- 
classical states. This approach or viewpoint allows us to 
give a rather detailed analysis of oscillations ||l^,0 in 
the PND, and to see precisely which kinds of oscillations 
stem from nonclassicality and which do not. As for the 
fact that we are able to obtain necessary and sufficient 
conditions for classicality, the principal tool we employ 
is the solution to the Stieltjes moment problem |^ in 
the classical theory of moments - when can a given se- 
quence of moments arise from a well-defined probability 
distribution? The Stieltjes problem deals with the case 
when the probability density is defined over the half line 
[0, c»). In our context, this corresponds to the range of 
the intensity variable of light. 

The material of the paper is arranged as follows. Sec- 
tion II collects basic definitions related to the diagonal 
coherent state representation of a single mode field, and 
discusses the classical - nonclassical divide from several 
angles. The final one chosen is that suited to phase insen- 
sitive measurements. Both the probabilities appearing 
in the PND, and the factorial moments of the PND, are 
recognised as moments of two auxiliary distributions over 
[0, oo), one with a clear physical meaning and the other 
a formal one. In Section III we give the basic method 
for obtaining local classicality conditions on the PND 
and obtain explicitly a minimal three term condition as 
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an example. This condition involves the photon number 
probabilities p„ for three successive values of n, in con- 
tradistinction to the condition in terms of the Mandel 
Q parameter which involves three successive (the zeroth 
moment which equals unity in every state, and the first 
and second) moments of p„ . The keen reader will notice 
indications of an interesting duality already at this ini- 
tial stage. We then put together a series of remarks and 
preliminary results based on these minimal conditions, 
which give a good orientation towards understanding the 
features of the PND from both local and global points of 
view. In particular we are able to discriminate between 
oscillations in the PND which genuinely reflect nonclas- 
sicality and others which do not. 

Section IV develops the complete necessary and suffi- 
cient conditions for classicality of the PND, by mapping 
the present problem on to the Stieltjes moment problem 
which is very well known and extensively studied. All the 
conditions obtained here are local in the sense mentioned 
above. From the necessary and sufficient conditions we 
then extract the conditions for classicality which go one 
step beyond the minimal ones discussed in Section II; this 
helps us sharpen some of the conclusions reached earlier. 

Section V gives an alternative or dual approach based 
on the factorial moments of the PND. This again is cast 
into the form of a Stieltjes problem, and the necessary 
and sufficient conditions for classicality are similar to the 
previous ones in structure. Section VI compares and con- 
nects the two treatments, and Section VII contains con- 
cluding remarks. 



II. BASIC DEFINITIONS, PRELIMINARY 
REMARKS AND RESULTS 

We consider a single mode quantized radiation field 
with photon creation and annihilation operators a) , a 
obeying the standard commutation relation 



[a, a}] — aa) — a) a = 1. 



(2.1) 



A general (pure or mixed) state of the field is described 
by a hermitian nonnegative density operator p with unit 
trace. According to the diagonal coherent state repre- 
sentation theorem, p can be expanded as an integral over 
projections on to the coherent states: 



p= _^(z)|z)(z| 

' TT 



(2.2) 



the integration being over the entire complex plane. Here 
the coherent state \z) is the normalized right eigenstate 
of a with (generally complex) eigenvalue z, related to 
the number operator eigenstates (Fock states) |n) in the 
standard manner: 



^)-xpf-v)f:-^ 



cxpl --|zr + 2a' ) |0) , 



\n) = (at)"|0)/Vn! ; 
a\z) = z\z) , a)a\n) = n\n) 



(2.3) 



The weight function 4'{z) in eqn.(2.2) is real and nor- 
malised to unit integral: 



(Pz 



(z) = l 



(2.4) 



However it is in general not pointwise nonnegative p3[ , 
and can be quite a singular quantity, namely a member 
of a certain precisely defined class of distributions over 
the plane. 

Any (hermitian) observable can always be written as 
a function of and d in normal ordered form, F {a^ , a) 
say . Its expectation value in the state p is then given 
by 



(F {a\a))^i:,{pF {a\a)) 



—<t>{z)F{z\z) 

TT 



(2.5) 



If in particular F (d^,d) is phase invariant, then its ex- 
pectation value does not require all the "information" 
contained in (/)(2;), and a simpler angle averaged auxil- 
iary distribution 7'(/) suffices: 

F(d^e^",de-^") = F {a\a) =^ 

OO 

(F(dt,d))= j dIV{I)F{l^/\l^/^) , 





27r 





OO 

j dIV{I) = 1. 



(2.6) 



We can regard 'P{I) as a real normalised marginal radial 
distribution function obtained from the complete ipiz). 
Note that -F(d^,d) being phase invariant is the same as 
it being a function of the number operator d^d. 

An important set of phase invariant F's leads to the 
photon number distribution (PND) or photon number 
probabilities {pn} in the state p: 

F„(dt,d) =: e-''^^^ : = |n)(nK 
Pn = {Fn {a\a)) 



ri=0 



/ 



dIV{I)e-'l"/nl ,n = 0,1,2,... . (2.7) 
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(The colons denote normal ordering). These p„ are al- 
ways well-defined for any bonafide p and always obey the 
laws for a discrete probability distribution: 

oo 

Pn>0, ^p„ = l. (2.8) 

ri=0 

With the two quantities 0(2;) and P(/) in hand, we can 
set up a three-fold classification of states p of steadily 
increasing nonclassicality. Conventionally p is said to be 
classical if (/)(z) itself is a probability distribution, that is, 
pointwise nonnegative and nowhere more singular than 
a delta function; otherwise it is nonclassical. In the for- 
mer case it follows that V{I) also can be interpreted as 
a probability distribution for the intensity. In the lat- 
ter case, a further refinement is possible based on the 
properties of V{I) and we arrive at the following scheme: 

p classical <^ 4>{z),V{I) > ; 
p weakly nonclassical <^ <j){z) ^0 , V{I) > ; 
p strongly nonclassical 4^ 0(z) /^O , V{I) /tO ■ (2.9) 

This is an exhaustive and mutually exclusive classifica- 
tion. The special role or significance of the PND {pn} 
can now be expressed as follows: even if V{I) is not 
a well-defined probability distribution and we are in the 
strongly nonclassical regime, the quantity V {I) always 
has well-defined finite moments n\p„ for all n. For ease 
in the following we introduce a special symbol for these 
moments, 

00 

g„ = nlp„ = J dlf{I) /" , 



V{I)=V{I)e-\ (2.10) 

the dependence on p being left implicit. On the other 
hand, the PND cannot discriminate between the classical 
and the weakly nonclassical cases in eqn.(2.9), in both of 
which V{I) > 0. Indeed, given a PND {p„} leading to 
a pointwise nonnegative V{I) (see below), if no phase 
sensitive quantities are to be measured we may take 

H^)^Vi\z\'), (2.11) 

and check that via eqn.(2.2) we get a physically accept- 
able density operator p. 

An interesting illustration of the weakly nonclassical 
case of eqn.(2.9) is the one-parameter family of states 
that results when a coherent state |zo) evolves through a 
Kerr medium for a time interval t. Indeed, for suitable t 
the state that results is the Yurke-Stoler state |^ 

\4,)^^{\zo)±z\-zo)). (2.12) 

This being the superposition of two coherent states has 
a (j){z) more singular than a tempered distribution. Nev- 
ertheless when angle averaged it leads to the same V{I) 



(and hence the same PND) as for the single coherent state 
|zo) : V{I) = 5{I - |zoP). The Hamiltonian of the Kerr 
medium being a function of a^a leaves the diagonal ma- 
trix elements Pn — {n\p\n) unaffected, but changes only 
the phases of {m\p\n) for m ^ n, and 'P{I) depends only 
on the diagonal elements of p and not on the off-diagonal 
elements. 

In the remainder of this work we shall be exclusively 
concerned with the PND {pn}, or more generally with 
expectation values of phase insensitive observables alone. 
For these purposes we shall combincdly refer to the classi- 
cal and weakly nonclassical cases of eqn.(2.9) as classical 
{V{I) > 0); and to the strongly nonclassical in eqn.(2.9) 
as nonclassical {V{I) ^0)- That is, for phase insensitive 
observables, we define: 

p classical <^ V{I) > ; 
p nonclassical P(/) ^0 . (2.13) 

Our principal aim now will be to develop necessary and 
sufficient conditions on the PND {p„} for p to be classi- 
cal; we may then say for brevity that we have a classical 
PND. 

The recent experiment of Munroe et al |l^ best illus- 
trates our considerations based on phase-insensitive mea- 
surements on the one hand and our classical-nonclassical 
divide (2.13) based on 7^(/) on the other. They show 
that the phase averaged quadrature amplitude distribu- 
tion P(^) measured by their optical homodyne detection 
apparatus is invertibly related to the PND. We shall show 
presently that the relationship (2.7) between P{I) and 
the PND is invertible. It thus follows that their proba- 
bility P(^) and our quasiprobability 1^(1) are invertibly 
related: their experiment extracts no more, and no less, 
information than contained in V{I). In other words, the 
classification (2.13), rather than (2.9), is the one relevant 
for such experiments. 

The distribution character of (j){z) passes over to a cor- 
responding property for V{I) - namely, it too is a member 
of a precisely defined class of distributions over [0,oo). 
While it is clear that the sequence {pn} cannot possibly 
capture all the information contained in ^(z) as all the 
off-diagonal matrix elements {m\p\n) , m ^ n, are ig- 
nored, we can easily show that the PND {p„} and the 
distribution P(/) determine each other uniquely. To re- 
cover V{I) from {pn\ we define a generating function 
A(A'), < X < 00, to represent the latter: 

00 

K{K) = Y,i-KrPn/n\ 

00 

= J dIV{I) e-fJoi2VlK) . (2.14) 


There is a great deal of freedom in the way we set up 
A(_Rr); we have chosen it so that, among other things, 
on the basis of eqn.(2.8) it is entire analytic in K. One 



3 



can then invert eqn.(2.14) by using the Fourier - Bessel 
integral theorem to get 



oo 

V{I) = J dK A{K)Jo{2VkI) . 



(2.15) 



Several signatures of nonclassicality of the PND are 
well known. The most familiar is the Mandel Q- 
parameter criterion which distinguishes (in a global 
sense) between super and subpoissonian PND's: 

Q = {{Anf - {n))/{n) 

= {{n') {nf {n))/{n) , 



oo „ 

^np„= / dIV{I)I , 

oo °? 

^nV. = / dIV{I)I{I+l) 

n 



classical PND Q > , superpoissonian statistics , 
Q < nonclassical PND, subpoissonian statistics . 

(2.16) 

However one can see that there are (uncountably many) 
states p for which Q is undefined, on account of either 
(n^) or both {n) and (n^) being divergent. This can 
happen because Q involves p„ for all n, and the rate of 
decrease of Pn as n oo may be insufficient for conver- 
gence of either (n) or (n^). A similar situation obtains 
with the factorial moments of the PND, namely 

7„, = (at"a™) 

En! 
in — m ! 

n—m ^ ' 

oo 

= Jdiril)!"', m = 0,l,,2,... . (2.17) 



In terms of these, the Mandel parameter is 

Q = (72 - 7?) /71 ■ (2.18) 

A well-known feature of the factorial moments is the fol- 
lowing iQ: 

1 /2 

Classical PND ^ 7,„'7m < -f,n' +ni < (72m' 72m) 

m', TO = 0,1, 2,... . (2.19) 

However, as with Q, 7^ also involves p„ for all but a 
finite number of values of n, and is undefined for the vast 
majority of states p. 

This discussion shows that in general, in contrast to 
V{I) = V{I) e^^, the moments 7™ of V{I) may not all 



exist. States p for which all 7^, or even all factorial 
moments upto some fixed maximum order, are finite are 
naturally quite severely restricted. 

One may enquire about the properties of the ordinary 
moments of the PND, namely 

Sm = {{a^d)'^) 

00 

00 



dIV{I) e-^^n™/"/ 



00 

= j dIV{I) e-' (^I ^ 



TO = 0,1,2, 



(2.20) 



While they may seem to behave qualitatively like the 7m, 
they are however not expressible as the moments of some 
function of / simply related to 7^(1) (as the p„ and the 
7m are, eqns.(2.7,10,17)). Therefore we do not deal with 
them hereafter. 



III. LOCAL CLASSICALITY CONDITIONS ON p„ 

The above discussion motivates the search for local con- 
sequences of classicality on the PND, that is, inequalities 
involving only a finite number of the p„'s which are al- 
ways well defined and which must be obeyed if the PND 
is classical. Based on eqn.(2.7), the general form of such 
conditions can be surveyed as follows. Consider a finite 
degree polynomial f{x) in a real nonnegative variable x, 
with the property of itself being real nonnegative: 



N 
n=Q 



> 0, for < X < c» 



(3.1) 



The nontrivial case here is when some coefficients c„ are 
negative; for example f{x) could be the square of a poly- 
nomial with coefficients of both signs. Then we find in 
view of (2.10): 



°? JV 
V{I) > ^ / dIV{I)e-'f{I) = ^ c„ g„ > 



. (3.2) 



One thus finds that for a classical PND, every polynomial 
obeying eqn.(3.1) leads to one inequality (3.2); so these 
are necessary local conditions for classicality. 

One can easily check that there are no such local condi- 
tions involving only two consecutive p„'s. The simplest or 
minimal local conditions involve three consecutive PnS, 
and are obtained by choosing for f{x) an integral power 
of X times a quadratic in the form of a perfect square: 
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Vil) > , fix) = - ax)2, a real =^ 

Qn-i ~ 2aqn + a^Qn+i > , all real a 

ql < qn-iqn+i , n = 1, 2, . . . . (3.3) 

Written in term of the true probabilities p„ this reads: 

1 



PND classical =4> p„ < (^1 + ~ ) Pn-iPn+i, 
n= 1,2,... 



(3.4) 



We shall call this sequence of local conditions the minimal 
or three-term classicality conditions on the PND. 

We now list a series of direct consequences of the above 
discussions and of the minimal local conditions (3.3) for 
a PND to be classical. 

(a) For a classical PND, the (nonnegative) distribution 
7^(7) is either concentrated at / = 0, V{I) — 5{I), 
or else it has positive definite weight for some values 
or ranges of / strictly greater than zero. The former 
corresponds to the vacuum state: 

V{I)=6{I)^Pn = qn=5^fi, 

P=|0)(0|. (3.5) 

In the latter situation we have 

V{I) > 0, ni) ^ 6(1) ^Pn>0, 



n== 0,1,2,. 



(3.6) 



Both (3.5) and (3.6) are immediate consequences 
of (2.7) or, equivalently, (2.10). Therefore we con- 
clude: (i) if po = 0, the PND is definitely nonclassi- 
cal; (ii) if the PND is classical and the state is not 
the vacuum, every p„ is strictly positive; equally 
well a non- vacuum classical state cannot be orthog- 
onal to any Fock state; (iii) conversely, for a nonva- 
cuum state, the vanishing of any one (or more) of 
the Pn's implies nonclassicality. 

These conclusions can also be obtained recursively 
from the minimal conditions (3.3): if Qno = for 
some uq > 1, repeated use of these inequalities 
leads to g„ = for all n > 1, the state being as- 
sumed to be classical. 

Hereafter we omit the vacuum state from consider- 
ation. 

(b) As a corollary to the above, we see that if p is any 
state and p' is obtained by adding some number m 
of photons to p, 



m > 1 



(3.7) 



then p' is always nonclassical. The normalisation 
constant M needs to be finite for (3.7) to make sense 
as a state. A direct calculation gives: 



M = tr(a™aT™p) 



in -f m)\ 



Pn = 



1 — Pn ; 

n=Q 

n < m — I 



(3.8) 



where p determines the PND {pn} and p' deter- 
mines {p'n}- If the series for J\f converges, then 
p' is a well defined state; and then the vanishing of 
p'qjP'i, ■ ■ ■ ,p'm~i establishes its nonclassicality. This 
proves that all photon added states are nonclassical. 
Nonclassicality of photon added coherent states jlj] 
and photon added thermal states [|6|p9|] (A/" is finite 
in both cases) have already been studied in great 
detail 

(c) For any classical PND, since go = Po > 0, we see 
that qo^V{I) is a nonnegative function normalised 
to unit integral. Therefore it can be treated math- 
ematically as though it were a probability distri- 
bution over [0, oo). However its physical interpre- 
tation is quite different from that of V{I) given 
earlier. 

(d) Since {qn} is a geometric sequence for a Poissonian 
distribution, we see that the minimal classicality 
conditions (3.3) are saturated in this case for ev- 
ery n. Thus we can interpret these conditions as 
requiring that the sequence {pn} be locally Poisso- 
nian or superpoissonian at each n. For a classical 
PND each q„ (respectively p„) has a lower bound 
determined by the values of the two previous q's 
(respectively p's): 

qn > ql^i/qn-2 , Pn > (^^ ^ ~ ^ Pn-l/Pn-2 , 

n>2 . (3.9) 

If even at one value of n > 1 we have 

ql > qn-iqn+i , (3.10) 

so that the PND is locally subpoissonian at that 
value of n, the state is definitely nonclassical. 

(e) The inequalities (3.3) also imply that for a classical 
PND we can never have g„ > qn+i , qn-i for any 
n > 1. That is, for classical states local maxima 
in the sequence {qn} o,re ruled out. This implies 
in particular that no oscillations in {qn} are possi- 
ble for such states. Thus we arrive at the impor- 
tant conclusion: oscillation in {qn} is a sure sign 
of nonclassicality. 

(f) On the other hand, the minimal conditions (3.3) 
permit a local minimum in {qn} but no more than 
one. This is because if we have minima at ni and 
again at n2 > ni + 2, there would necessarily be a 
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local maximum in between, which is disallowed for 
classical states. 
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FIG. 1. "Classical" oscillation in the photon number dis- 
tribution, for an incoherent mixture of coherent states with 
density matrix p — Ai|qi) (qi 1+ A2|a2){a2|+ A3|a3){Q:3|-|- 
A4|a4)(a4|+ A5|a5)(a5| with |ai|^ = 10, |a2p = 30, 
laal^ = 60, |a4p = 90, [asl^ = 130 and Ai = A2 = 0.25, 
A3 = 0.2, A4 = 0.18 As = 0.12. The symbols p{n), q{n) stand, 
respectively, for Pn,qn of the text. It may be noted that q„ 
exhibits no oscillations for this classical state. 



(g) We can now combine (e) and (f) above to state the 
following. Ever since the important work of Schle- 
ich and Wheeler on interference in phase space , 
the statement that oscillation in {pn} is a signature 
of nonclassicality has almost become a folklore. In- 
deed, oscillations in {p„} are known as nonclassical 
oscillations The striking virtue of this char- 

acterization is that it is local in n, in the spirit of 
our present approach, as against other characteri- 
zations based on the (factorial) moments of {pn}- 
We have shown in Fig.l the sequence {pn} for a 
suitably chosen incoherent superposition of coher- 
ent states. This state is classical by construction, 
yet it exhibits oscillations in {pn}, showing that the 
above characterization needs quantification of some 
sort while retaining the attractive feature of being 
local in n. Our minimal local conditions (3.3) can 
be viewed as a quantification of this type. They 
limit the extent to which oscillations can occur in 
{Pn} if the state is classical. These limits are ob- 
tained by translating properties (e), (f) of {g„} 
above into corresponding properties of {pn}- For 
classical {qn}, (e) and (f) allow only four generic 
patterns of behaviour: (i) q„ nondecreasing, pro- 
vided a constant phase (if any) is followed, not pre- 
ceded, by an increasing phase (if any); (ii) q.^ non 
increasing, provided a constant phase (if any) is 
preceded, not followed, by a decreasing phase (if 
any); (iii) g„ constant (actually subsumed in the 



two previous possibilities); (iv) q„ passing through 
one local minimum, monotonic decreasing (increas- 
ing) before (after) the minimum. The monotonic 
increasing case is part of (i) , the monotonic decreas- 
ing case is part of (ii). Well known classical states 
constitute examples of these generic behaviours of 
classical {qn}, as shown in Fig. 2. 




n 

FIG. 2. Generic behaviours of the sequence {qn} for dif- 
ferent classical states. First three are coherent states with 
different displacements and the fourth one is a thermal state. 

It turns out that constant phases in {qn} in cases (i) 
and (ii) above, allowed by the minimal local conditions 
(3.3), but not shown in Fig. 2, are indeed forbidden by 
the higher order local conditions as will be shown in Sec- 
tion IV. 

All these restrictions on classical {qn} translate into al- 
lowed behaviours for classical {pn}- First, local maxima, 
Pn > Pn±i, are permitted, provided 



Pn 



Pn 



< 1 



1 



Pn~l Pn+1 



(3.11) 



Next, again because of the factor (l + ^) on the right 
hand side in the minimal classicality condition (3.4), 
some amount of oscillation in {pn} is allowed. This is 
just the kind of oscillation seen in Fig.l. Note that the 
period of oscillation (difference of n values at two suc- 
cessive maxima in {pn}) in Fig.l is substantially greater 
than two as against the period two oscillations occurring 
in the PND of a squeezed or cat state. For period two 
oscillations the conditions (3.4) place substantial restric- 
tions on the amplitude. (Indeed the next higher order 
local conditions to be derived in Section IV make these 
restrictions even more stringent). Further, as the factor 
(1 -I- i) approaches unity as n becomes large, the restric- 
tions on the amplitude of a classically allowed period two 
oscillation are stronger at higher values of n. To end the 
discussion in this paragraph we note that the period two 
classical oscillation in {pn}, allowed by our lowest order 
local condition (3.3,4), indeed survives all higher order 
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local conditions, as can be seen by the example shown in 
Fig.3. 




n 

FIG. 3. Limits to PND oscillations from classicality con- 
ditions- an example of incoherent mixture of three coher- 
ent states: p = Ai|ai)(ai| + A2ja2)(a2| + A3|a3)(a3| with 
\aif = 0.0, |q2|^ = 1.56, |a3p = 5.42, Ai = 0.06, A2 = 0.305 
and A3 = 0.635. Period two oscillations though restricted in 
amplitude can be present for a classical state. 

An instructive example illustrating many of the points 
made above is provided by the case of a pure state ob- 
tained as the superposition of two "opposite" coherent 
states: 



AA = 2 (l + e-^l^ol'cos^) . (3.12) 



Here 9 is the relative phase (in the Pancharatnam 
sense pT[) between the two components of the super- 
position. As special cases we get the Yurke-Stoler states 
(2.12) when 6 = ±7r/2, and the cat states when 6 — Q,7t. 
For the PND we have 



p„ = |(n|V'(zo,0))P 

_.-UoP l^ol'" (1 + (-!)" cos( 



1 + e~2|^"l' cos 



(3.13) 



so that 



/ 2 

9ri-l9n+l/9„ 



{(1- 
fniO) 



-l)"cos6l)/(l + (-l)"cos6')}^ 

(3.14) 



Clearly, /„(0) < 1 for even n if -7r/2 < 9 < 7r/2; and 
fn{9) < 1 for odd n if 7r/2 < 61 < 37r/2. Thus the state 
\'4'{zq, 9)) violates the minimal local condition (3.3) and is 
nonclassical for all ^ ±7r/2. For 9 = ±7r/2 these condi- 
tions are saturated at every n, and we have a Poissonian 

{Pn}- 

It is clear that nonclassicality of the Yurke-Stoler state 
(2.12) can be exhibited only through phase sensitive 



considerations, since the PND is poissonian. It is also 
well known |1^ that all pure states other than coher- 
ent states are nonclassical at least at the phase sensitive 
level. What is really interesting is the fact that for every 
9 ^ ±7r/2, the nonclassicality of |V-'(-zo, ^)) is coded in the 
phase insensitive PND, and that our minimal conditions 
(3.3,4) capture this nonclassicality! 

We conclude this Section with one more application of 
the local condition. Recent years have witnessed a re- 
markable progress in quantum state reconstruction using 
techniques of optical homodyne tomography [ p2[ . As a 
consequence it is now possible to 'map out' the Wigner 
distribution of a state using the inverse Radon transform, 
or reconstruct the density matrix in the Fock basis us- 
ing a set of pattern functions. Schiller et al report 
such a reconstructed density matrix pm,n for m, n < 6. 
The reported values of g„ for n = to 6 are 0.44, 0.07, 
0.26, 0.30, 1.44, 3.60 and 28.80. The local conditions 
are clearly violated: qiq^ < g|, q^q^ < q\. Hence the 
state is nonclassical. The state was, of course, known to 
be quadrature squeezed, and hence nonclassical. What 
is interesting about the present illustration is the fact 
that the nonclassicality of the Schiller et al state sur- 
vives phase averaging, and that we are able to arrive at 
this definitive conclusion with just a few values of the 
diagonal elements of the density matrix! 



IV. NECESSARY AND SUFFICIENT 
CONDITIONS FOR CLASSICALITY OF THE 
PND - THE STIELTJES MOMENT PROBLEM 

We have seen in the previous Section that pointwise 
nonnegativity of 7^(/), hence of 7-'(/), leads to the mini- 
mal local conditions (3.3,4) on the sequences {g„}, {pn}- 
As already remarked, these are necessary conditions for 
classicality. Central to this derivation was appreciation of 
the fact that the g„ are moments of V (1) , eqn.(2.10). We 
have also seen that if the state is classical, then q^^V^) is 
well-defined and mathematically interpretable as a prob- 
ability distribution over [0, 00). In this Section we exhibit 
the necessary and sufficient conditions on the PND {qn} 
in order that the state p be classical. We then use these 
conditions to examine the extension of the inequalities 
(2.. 23, 24) to the five-term case. 

The reconstruction of a probability distribution from 
its moment sequence constitutes the classical moment 
problem on which there exists an enormous amount of 
literature [T^ ] . When the probability distribution is over 
[0,00), one calls it the Stieltjes moment problem. The 
Hamburger moment problem corresponds to the case 
where the probability distribution is over (—00,00); and 
the Hausdorff or the little moment problem corresponds 
to the finite interval [0, 1]. Since the argument of 7^(1) is 
nonnegative, our problem of deriving necessary and suf- 
ficient conditions on the moment sequence {qn\ in order 
that qQ^Vi^I) is a true probability distribution is then a 
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Stieltjes moment problem. In this Section we can take 
advantage of the fact that these moments {q-n} of V{I) 
are always well-defined and finite, whatever the state p 
may be. 

The solution of this classical problem is well known. 
To exhibit it, construct from the sequence {g„} two sym- 
metric matrices 

L^N)^l{N) dimension (A^ 1) defined 

by 



^mn Hm+n i 

I^iZI = 9m+n+l , TO, n = 0, 1, 2, . . . , N. 



(4.1) 



That is, 



l(N) 



( 10 



qi 
q2 



92 

qs 



qN \ 
qN+i 



\ qN qN+i qN+2 ■ ■ ■ qiN J 
( 



qi 
q2 



92 
93 



93 
94 



9JV+1 \ 

9JV+2 



(4.2) 



\ 9^+1 qN+2 qN+Z ■ ■ ■ 92JV+1 / 



Thus, L^^) arises from 

2^(^+1) by deletion of the first 
column and the last row in the latter. Notice also that 
the diagonal elements of L^^'^ are qn for even n < 2N, 
while those of Z^^-' are g„ for odd n < 2N + 1. 

Now we quote the fundamental classical theorem . 
Theorem 1: 

The necessary and sufficient condition on the PND 
sequence {qn} in order that the associated distribution 
qQ^V{I) be a true probability distribution over [0,oo) is 
that the above matrices be nonnegative for all iV > 0: 

Classical PND ^ P(/) > ^ P(/) > 

^LW,iW >0, 7V = 0,1,2,.... (4.3) 

Proof: 

We make use of the following two facts. 

(i) A distribution V{I) over [0, oo) is pointwise nonneg- 
ative if and only if it leads to a nonnegative expectation 
value for every polynomial /(/) which is itself pointwise 
nonnegative over [0,oo) p4[ : 

oo 

Vil) > for aU / ^ (f{I))p - / dIV{I)f{I) > 0, 



for aU polynomial /(/) > 0. (4.4) 

(ii) Any polynomial /(/) pointwise nonnegative over 
[0, c») can be written in terms of two perfect square poly- 
nomials as [|2| 

f{i) = {Mi)r+i{.f2{i)r ■ (4.5) 

To prove necessity, suppose that V{I) is pointwise non- 

N 

negative. Consider the polynomial /i(/) = ^^nl^, 

n=0 



where the c„ are arbitrary real coefficients. By (4.4) we 
have {[h{I)f)f > 0. That is, 

N 

m.n— 

N 

m,n— 
TV 

= E c™c„lW>o. (4.6) 



m.n— 



This means L(^) > for every iV > 0. Similarly, writing 

N 

/2(/) = J2 dnl"" and evaluating the expectation value 

of the nonnegative polynomial /(/2(/))^ we have from 
(4.4): 

N 
m,n— 
— ^ ^ dm dn Qrn + n + 1 



m,n— 
TV 



(4.7) 



m,Ti— 



This means L^^^ > for every iV > 0. Thus nonnega- 
tivity of P{I) implies that of every L^'^\L^^\ 
To prove sufficiency, assume 

Li^)^Ui^) >0for7V>0. 
Given any nonnegative polynomial /(/), writing it in the 
form (4.5) we find that {/{!)) f, > 0. This implies, by fact 
(i) stated at the beginning of the proof, that 'P(/) > 0. 
This completes the proof of the Theorem. 

Now each matrix L^^^ (respectively L^^^) arises from 
7^(^-1) (respectively L^^~^^) by addition of one extra 
row and column, leaving the rest intact. Therefore the 
result (4.3) translates into a series of determinant condi- 
tions: 



Dn = detL^^) , Dn = detL(^) : 
Classical PND ^V{I),V{I)>0 

^ Dn , Dn>0 , iV = 0,1,2, 



(4.8) 



Assume now that we have a classical PND. Then either 
Dn, Dn>0 for all N, or Dn, Dn > for N < k and 
Dn = Dn = for > fc. In the latter case the support 
of V{I) (the set of values of / at which ?-'(/) > 0) is a 
finite set of k points |2^]. In the former case it is an 
infinite set. This can be intuitively understood along 
the following hues. Suppose V{T) = S{I — Iq). Then 
qn - /o" so that L^„^ = I^I^ and = WIq- That 
IS, L(^) and L(^) are (essentially) projection matrices. 
Thus when the support of V{I) is a finite set of points, 
L^^^ (as also Z*^^)) is the sum of k projections. Since 
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/o = contributes a projection only to L^^^ but not to 
, one has the following refinement: if the support of 
V{I) consists of k points including the point / = 0, then 
Dn > 0{Dn > 0) if and only if TV < k{N <k-l). 

It is useful to remark that the first fact (4.4) we have 
used in the proof of the theorem is common for all the 
three types of moment problems. What changes from one 
moment problem to another is the second fact dealing 
with the decomposition of nonnegative polynomials into 
sums of square polynomials, thus enabling us to convert 
(4.4) into simple matrix conditions. For instance, for 
the Hamburger moment problem on (—00,00) we have, 
in place of (4.5), the statement that every polynomial 
nonnegative over (—00,00) can be written as the sum 
of two square polynomials. Thus in this case we have 
to deal only with the matrix Li^\ and the conditions 
L(^) > for aU TV > are both necessary and sufficient 
for positivity of the underlying Hamburger distribution. 

It is immediate to relate the minimal local condi- 
tions (3.3,4) to the above theorem. Nonnegativity of 
L<-N)^L'-N) demands as a necessary condition nonnega- 
tivity of every 2x2 determinant obtained from 2x2 
blocks along the principal diagonals of L^^^ and U^K 
This is precisely the condition (3.3). It is also clear why 
these minimal conditions are only necessary: positivity of 
the diagonal 2x2 blocks of Li^\ L^-'^^ does not capture 
in its entirety the positivity of L^^^ and Li^\ 

We now go beyond the minimal conditions (3.3,4) and 
derive the next hierarchy of local conditions for a classical 
PND. Given {(?„} we define 



1^1 



1,2, 



(4.9) 



This definition is legitimate since in any case for a classi- 
cal PND each g„ > 0. Then the minimal conditions (3.3) 
simply read: 



Classical PND ^ a;„ > 1 , n = 1, 2 . 



(4.10) 



As one may anticipate, the next order local conditions 
which we now derive, and which are still only necessary 
conditions for classicality, involve five consecutive gn's, 
or cquivalcntly three consecutive x„'s. A necessary con- 
dition for the nonneg ativity of L(^),L(^) is that their 
diagonal 3x3 blocks be nonnegative. That is, 



after elementary algebra the third reduces to a condition 
expressible in terms of three consecutive x„'s: 



n+1 



1) > 



1 



n = 2, 3, 



(4.13) 



These are our next to minimal or second order local con- 
ditions on the PND for classicality. 

To conclude this Section, we present an interesting 
implication of these conditions. As already noted, if 
{pn} is Poissonian then {qn} is a geometric sequence, 
so x„ = 1 identically and (3.3) are saturated. We now 
ask whether it is possible to have a classical state for 
which qn-iln+i = Qn for some values of n > 1, whereas 
Qn-ign+i > Qn ^0^' Other n's. Such classical states, if they 
exist, can be said to be locally Poissonian at the former 
values of n, namely whenever a;„ = 1. 

Suppose a classical state is locally Poissonian at no, 
with — 1. Then two applications of (4.13), once with 
n = uq — 1 and then with n — iiq + 1, lead to Xno±i = 1. 
Continuing this process we end up with a:„ = 1 for all n > 
1. Thus a classical state cannot be only locally Poissonian 
: it is either locally Poissonian throughout, (.t„ = 1 for 
all n > 1) or locally superpoissonian throughout {xn > 1 
for all n). Translating this to {pn} we can now strengthen 
our minimal conditions (3.4) to say: 

For a classical PND, either 

Pn-iPn+i = —^Pn (Poissouian) , n > 1 
n+1 

or 

Pn-iPn+i > ^ , Pn (Supcrpoissonian) , n > 1. (4.14) 
n -f 1 

This is a refinement of the minimal or first order neces- 
sary local classicality conditions (3.4) in the light of the 
second order conditions (4.13). 

It is now clear that the constant phases allowed in pat- 
terns (i) and (ii) under (g) preceding eqn.(2.31) are for- 
bidden by the above refined condition (4.14). This ren- 
ders Fig. 2 generic and exhaustive. 



ALTERNATIVE FORMULATION VIA 
FACTORIAL MOMENTS 



9n-2 Qn-l In 
Ai= I Qn-l qn Qn+l | > , n = 2, 3, . . . . (4.11) 

Qn Qn+1 Qn+2 

If we bring in Xn-i, Xn, Xn+i here and eliminate qn±2 h^ 
favour of the other variables, this reads: 

I I Ql-lXn-l Qn-lQn qI 

An^ Qn-lQn Qn QnQn+l | > 

\ qI QnQn+l Ql+iXn+l 



n = 2,3,. 



(4.12) 



This entails three conditions: q^_iXn-i > 0, Xn-i — 1 > 
0, detA„ > 0. The first two are already obeyed; and 



In this Section we present an approach to classicality 
of a PND based on the normal ordered moments of p 
or factorial moments of {pn}, namely the sequence {7m} 
defined in cqn.(2.17). This approach will be along the 
lines of earlier work of Agarwal and Tara however 
our conditions for classicality will be both necessary and 
sufficient. As is amply clear from the discussion in Sec- 
tion II, though, we must bear in mind that for the vast 
majority of states p and PND's {pn}, the 7™ are not all 
defined. 

Suppose we have a PND {pn} whose factorial moments 
{7m} are all finite and known. (This certainly places very 
severe restrictions on the behaviour of Pn for large n). 
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Our problem is to find necessary and sufficient conditions 
on {7m} in order that the PND be classical. Now the 7™ 
are the moments of just as in the preceding Section 
the Qn were moments of V{I) — V{I)e~^ . At the risk 
of repetition, the finiteness of all 7™ places very strong 
restrictions on 

It is now clear that we have again a Stieltjes moment 
problem. We want the necessary and sufficient conditions 
on {7m} to ensure pointwise nonnegativity oiV{I). Form 
two matrices mW,MW of dimension using the 

7's: 



This leads to the following matrices M'^^^ and M*^^^: 



/ 70 

71 



71 

72 



72 

73 



\ 7A' 7^+1 lN+2 

/ 71 72 • ■ • 

72 73 • ■ ■ 

\ 7JV+1 lN+2 ■ ■ ■ 



IN \ 

12N ) 

7^+1 \ 
7JV+2 

72Af+l / 



(5.1) 



Both are real symmetric; A/^^^ equals M*^^+-'^) with first 

column and last row deleted; and M^^^ ^iVf^^^^ has 7m 

for m even (odd) along the diagonal. Thus, exactly as in 
Section IV, we have: 
Theorem 2: 

The necessary and sufficient conditions for the PND 
{p„} to be classical (given the existence of all factorial 
moments 7^) are 



> 0, M(^) > 0, TV = 0, 1, 2, 



(5.2) 



The proof is exactly parallel to that in Section IV, 
with ■P(/), L^^^, LS^^ replaced respectively by 

V{I), 7«, MW, M(~). 

This theorem completes the work initiated by Agarwal 
and Tara |^ , by improving their necessary conditions for 
classicality (they had only the condition M^^-* > 0) into 
the necessary and sufficient condition (5.2). Thus the 
constraints on the factorial moments arising from the re- 
quirements M*^^^ > are the same as in their work: 
with TV = 1 we have 7072 — 7i > 0, ie., 72 — 7i > 
which (see eqn.(2.18)) is the same as requiring the Man- 
del Q-parameter to be nonnegative. With TV = 2 we 
obtain the additional condition detM^^^ > 0, and so on 
as in H . However the constraints arising from M(^) > 
are new. For TV = we have 71 > 0; for TV = 1 we have 
71 73 — 72 ^ 0; and so on. 

To conclude this Section, let us re-examine the class of 
pure states |-0(zo,6') > defined in eqn.(2.32) within the 
present approach. All the 7m do exist in this case and 
are easily calculated. 



7m = (V'(zo,e)|at™a'"|V(zo,^^)) 

(^l + (-l)™e-2|^oPcos6') 



l^ol 



((1 



-2|zo|- 



COS! 



(5.3) 



TW(^) = A 



TW(^) = B 



1 


a 


1 


a . . . 


a 


1 


(7 


1 . . . 


1 


a 


1 


a . . . 


a 


1 


a 


1 ... 


1 


a 


1 


a . . . 


a 


1 


a 


1 ... 



A 



B, 



^-diag(l,|zo|M^olM^o|',. 
B = diag(|zo|,|zo|Mzo|^...) 

\ - e-^l'^ol' cos( 



-2|zo|- 



(5.4) 



It is clear that both Af and TVf(^) are matrices of rank 
2. It is further clear from their structures that TW^^^ > 
if and only if cr < 1, while TVf(^) > if and only if cr > 1. 
We thus get again the same results as in Section II: for 
7^ ±7r/2, the state \'il){zQ, 9)) has a nonclassical PND. 

We wish to remark once again that the fact that the 
conditions M^-^^ > in (5.2) are needed over and above 
Tkf(^) > to form a set of necessary and sufficient 
conditions for classicality is ultimately due to the fact 
that our moment problem is a Stieltjes problem: had it 
been a Hamburger problem, the Agarwal- Tara conditions 
M > would have been both necessary and sufficient! 



VI. CONNECTION BETWEEN THE TWO 
APPROACHES 

We have presented two approaches to the problem of 
phase - insensitive nonclassicality of a quantum state p: 
one based on {(?„}, the moment sequence of V{I); the 
other dual approach based on {7m}, the moment se- 
quence of V{I) when defined. In each case we obtained 
necessary and sufficient conditions for a PND to be clas- 
sical, exploiting the fact that the underlying problem was 
a Stieltjes moment problem. In this Section we bring out 
explicitly the connection between these dual approaches 
and establish their equivalence, in the case (naturally) 
when all 7m are defined. 

The fact that 7^(1) — V{I)e^^ suggests the use of 
Laplace transforms Let $(s),$(s) be the Laplace 

transforms of V{T),V{T) respectively: 



00 

$(s) = J dl 'Pil) e" 



$(s) ^ J dl V{I) e 


= $(s + 1) . 



si 



-si 



(6.1) 



10 



Wc now exploit the fact that the moments of a distribu- 
tion are simply related to the derivatives of its Laplace 
transform at the origin: 

oo 

7n = ldini)r = {-ir^l^,; 



oo 

,„ = /d7P(7)e-7" = (-ir^U 



(6.2) 



Making two Taylor series expansions of $(s), once about 
s = and then about s = 1, equating the two expansions 
and using eqn.(6.2), we have: 

oo oo 

= E(-l)'«^(^ - 1)V^! • (6-3) 



fe=0 



Equating derivatives of both sides first at s = and then 
at s = 1 gives the pair of relations 



In = ^Qn+k/kl , n > ; 

fe=0 
oo 

qn = , n > . 



(6.4) 



fe=0 



Writing the two sequences {(/„}, {7n} a-s two column vec- 
tors Q and T respectively we have the infinite matrix 
equations 

r = SQ, 

1 1/1! 1/2! 1/3! ... 
1 1/1! 1/2! ... 
1 1/1! ... 



1 -1/1! 1/2! -1/3! ... 
1 -1/1! 1/2! ... 
1 -1/1! ... 



(6.5) 



We have displayed S, to exhibit the fact that these 
matrices are unchanged if the first row and first column 
are deleted. We shall have occasion to return to this 
important feature. 

While Q and F considered as (infinite dimensional) col- 
umn vectors are connected by the matrix S, the corre- 
sponding infinite dimensional matrices L = L^°°^ and 
M = M^°°^ are connected through a symmetric transfor- 
mation using 5^/^. We have: 



(5V2)., = 2^--V(fc-i)!, 



(So both matrices arc upper triangular: these matrix ele- 
ments vanish for k < j). That S^^^ and 5*^^/^ so defined 
are indeed inverses of one another follows simply from the 
familiar properties of binomial coefficients. Using these 
same properties it may be verified that 



L = S-^'^M{S-^'Y 



(6.7) 



(6.6) 



This proves that M > if and only if L > 0. 

Let Q be the moment sequence derived from Q by sim- 
ply dropping go- Clearly, L = L^°°^ is in the same rela- 
tion to Q as 1/ is to Q. Similarly if we form T from V 
by simply dropping 70, then M = M(°°) is in the same 
relation to F as M is to F. Now recalling the fact that 
5, S'^/^, 5"^/^ have the interesting property of each 
being unchanged upon deletion of the first row and first 
column we conclude: 

f = 5g , Q = S-^f ; 

M = S^'^L{S^'Y , L = S-^/^M{S-^/Y ■ (6.8) 

This proves that M > if and only if L > 0. Com- 
bined with the statement after eqn.(6.7), we have thus 
established the equivalence of the two approaches. They 
are dual to one another, and in particular Theorem 1 is 
equivalent to Theorem 2. 

The above analysis shows that if wc drop the first k 
terms from the moment sequence of a bonafide Stieltjes 
probability distribution (in the semi- infinite interval), the 
result is again a bonafide Stieltjes moment sequence (of 
some other valid probability distribution). This is a dis- 
tinguishing feature of the Stieltjes moment problem. It 
is not difficult to see that the corresponding statement 
is true for the Hamburger moment problem (on the en- 
tire real line) only if an even number of initial terms are 
dropped from the moment sequence. 



VII. CONCLUDING REMARKS 

Every probability distribution {pn} over the nonnega- 
tive integers, admissible according to the laws of classical 
probability theory, can appear as the PND of some quan- 
tum state p of the single mode radiation field. Wc have 
given explicit necessary and sufficient conditions - in two 
equivalent or dual forms - for a PND to be the result of 
a classical state in the sense appropriate to quantum op- 
tics. From the perspective of classical probability theory, 
the significant points are these: that we arc able to view 
the quantities g„ = n!p„ as the moments, in the Stieltjes 
sense, of an auxiliary function V{I), and that the prop- 
erties of P{I) or equally well of Vil) determine whether 
the state is classical or not. The route from {pn} to V{I) 
is via the pair of eqns.(2.14,15). 

We have emphasized the following points in our work: 
(i) that local classicality conditions on the PND are nat- 
urally available for all states p; and (ii) all statements 
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based on the factorial moments make sense only for a 
quite limited set of states p, for which the probabilities 
Pn go to zero "sufficiently fast" as n — > cx). It seems to 
us that some of these points have not been given in the 
past as much attention as they deserve. Above this, we 
stress the completeness of our results, and the more care- 
ful treatment of oscillations in the PND as a signature 
of nonclassicality, which becomes possible with our local 
approach. 
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